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ABSTRACT 

We present updated results of a high-resolution, magnitude limited (K^ < 
12 mag) imaging survey of nearby low-mass M6.0-M7.5 field stars. The obser- 
vations were carried out using adaptive optics at the Gemini North, VLT, Keck 
II, and Subaru telescopes. Our sample of 36 stars consists predominantly of 
nearby (< 30 pc) field stars, 5 of which we have resolved as binaries. Two of 
the binary systems, 2MASSI J0429184-312356 and 2MASSI J1847034+552243, 
are presented here for the first time. All 5 discovered binary systems have 
separations between (/.'OS — 0'.'53 (2-9AU) with similar mass ratios (q>0.8, 
AKg < Imag). This result supports the hypothesis that wide (a> 20 AU) very 
low-mass {Mtot < 0.19 M©) binary systems are rare. The projected semimajor 
axis distribution of these systems peak at ~ 5AU and we report a sensitivity- 
corrected binary fraction of 9^3% for stars with primaries of spectral type M6.0- 
M7.5 with separations > 3 AU and mass ratios q> 0.6. Within these instrumen- 
tal sensitivities, these results support the overall trend that both the semimajor 
axis distribution and binary fraction are a function of the mass of the primary 
star and decrease with decreasing primary mass. These observations provide 
important constraints for low-mass binary star formation theories. 

Subject headings: instrumentation: adaptive optics — binaries: general — stars: low 
mass — stars: individual () 
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1. Introduction 

One of the main motivations for measuring the binary fraction of stars is to better under- 
stand the process of star formation itself. After all, stars like our own Sun are preferentially 
produced in multiple systems (Duquennoy & Mayor 1991). The classic stellar formation 
mechanism of molecular cloud core collapse and fragmentation, however, has a hard time 
explaining the tightest systems. While this mechanism can explain wide binary systems 
(semimajor axis > 10 AU), it has some difficulties explaining the formation of tight systems 
(Bate, Bonnell, & Bromm 2002). Additionally, the multiplicities of the lowest mass stars 
and brown dwarfs appear to be statistically different from those of more massive systems 
(Close et al. (2003) and references within). These differences, if proven to be real, provide 
important clues and constraints to theoretical stellar formation models. 

The continuously improving statistics of binary stars brings clarity to the paradigm that 
the binary fraction and semimajor axis distribution are functions of the central star mass. 
Surveys of G dwarfs estimate a multiplicity fraction of approximately 50% for separations 
of 3AU and greater (Duquennoy & Mayor 1991). Other surveys of similar sensitivity to 
systems wider than 3AU have found that early M dwarfs (M0-M4) have measured binary 
fractions of ~ 32% (Fischer & Marcy 1992) while late M/early L dwarfs (M8.0-L0.5) estimate 
fractions ~ 15% (Close et al. 2003). The trend appears to continue to the coolest objects 
- L dwarfs reporting ~ 10-15% (Bouy et al. 2003; Gizis et al. 2003) and T dwarfs at 
~ 10% (Burgasser et al. 2003). The same surveys infer semimajor- axis separations to also 
be a function of primary mass. While G and early M dwarfs (M0-M4) show broad separation 
peaks of ~ 30 AU, late M (> M8), L, and T dwarfs appear to have separations tightly peaked 
at ~ 4AU (Close et al. 2003). A similar result has been shown to apply to the sequence of 
members in the Pleiades cluster, from solar-type stars to brown dwarfs (Martin et al. 2003). 
Together, these results are providing both clues and empirical constraints on star formation 
mechanisms as well as potentially help calibrate the mass- age-luminosity relation for very 
low-mass (VLM) stars. 

In Siegler et al. (2003), hereafter referred to as Paper I, we presented results from the 
largest flux hmited (K^ < 12 mag) survey of nearby fleld M6.0-M7.5 dwarfs. The binary 
fraction of this narrowly deflned spectral type range had not been quantifled as those of 
stars slightly earlier (M0-M4; Fischer & Marcy 1992) and later (M8-L0.5; Close et al. 2003). 
Considering the differences in binary charactersitics as a function of mass as discussed above, 
would M6.0-M7.5 binaries have intermediary characteristics to their main sequence neighbors 
or would they resemble the ultracool dwarfs? 

Paper I's sample consisted of 30 stars and presented the discovery of three new binary 
systems using the University of Hawaii visitor AO system Hokupa'a (Graves et al. 1998) at 
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the Gemini North telescope. The discoveries followed characteristics of other VLM binary 
systems, namely relatively equal mass components (q>0.8) with projected separations less 
than 16 AU. In this paper we present our two latest binary discoveries from this spectral 
range, 2MASS1 J0429184-312356 and 2MASS1 J1847034+552243, hereafter referred to as 
2M 0429 and 2M 1847. These binaries were discovered with the VLT and the Subaru AO 
facilities, respectively. The total M6.0-M7.5 sample size is increased to 36 and we update 
the binary fraction results with those presented in Paper I. We present our observations 
and results in the following section and examine the systems' derived characteristics such as 
distance, age, temperature, spectral type, and mass in §3. In §4 we conclude by discussing 
the binary frequency and separation distribution of M6.0-M7.5 dwarfs. 



2. Observations and Results 
2.1. The Sample 

We selected a flux-limited sample of 36 objects consisting of M6.0-M7.5 dwarfs with 
Kg < 12 mag and J-K^ > 0.95 mag from mainly 2MASS stars listed in Cruz et al. (2003), 
Reid et al. (2002), and Gizis et al. (2000). Paper I discusses the first 30 observations and we 
report here the most recent 6. One of the 6 targets is a recently discovered high proper motion 
M dwarf - SO 025300.5+165258 (~ 3.6 pc away; Teegarden et al. 2003). We discuss this star 
further in §2.5. We also note that we observed at Subaru on 2003 July 10 (UT) another 
recently-discovered high proper motion M dwarf 2MASSI J1835379-F325954 (~ 5.7 pc away; 
Reid et al. 2003). While not part of this sample due to its later spectral type (M8), it was 
observed at high resolution and found to have no q> 0.8 companions at separations > Of!l. 



2.2. The Telescopes and their AO Systems 

The 30 targets from Paper I were all observed at the Gemini North telescope. The 6 
targets presented here were observed at the Subaru and VLT Observatories. Due to its recent 
discovery and proximity, we conducted additional long integrations of SO 025300.5-1-165258 
with the Keck II telescope. Interestingly, the AO systems at these telescopes represent the 
three major wavefront sensor (WFS) technologies currently in use today. Gemini North, at 
the time of our observations, and Subaru, currently, use 36-element curvature WFSs, the 
VLT has an infra-red Shack-Hartman, and the Keck II utilizes a visible Shack-Hartman. This 
survey provided the opportunity to compare and contrast how different AO WFSs differ in 
their abilities to lock on faint targets. 
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As discussed in detail in Paper I, one of the challenges in utilizing AO is locating 
sufficiently bright guide stars near enough one's science objects to minimize uncertainty in 
the image quality introduced by isoplanicity. This is best achieved when using the target 
object itself as the AO guide star. This shifts the criterion of target selection from the 
availability of bright natural guide stars (R<15 mag) to the sensitivity of the respective 
telescope's AO system, in particular the WFS. This becomes quite important because the 
probability of locating a R=15mag star within 30"of one's science target is only about 
~ 15% (Fig. 3.10; Roddier 1999). The ability to guide on fainter stars also allows for both 
larger sample sizes and improved contrast ratios. 

We were able to observe the faintest of our targets (V~19.0- 19.5 mag, 1~15.5- 16 mag) 
only with the former Gemini North AO system Hokupa'a where we conducted the ma- 
jority of the observations. No other current AO system can lock onto such faint targets. 
Hokupa'a, decommisioned in 2003, was a curvature-based AO system which employed in its 
WFS red-sensitive, photon-counting avalanche photodiodes with effectively zero read-noise. 
Consequently, this type of sensor is ideally suited for guiding on intrinsically faint objects 
as long as they are relatively red (V-l~ 4 mag). The Keck 11 telescope with a more tradi- 
tional Shack-Hartman WFS allows for improved angular resolution with higher obtainable 
Strehl ratios but requires brighter targets (V~ 15 mag). We compare and contrast the per- 
formance of these two types of WFS technologies in Siegler, Close, & Freed (2002). At both 
Subaru and the VLT we were able to lock on our faint low mass targets with I< 15.2 mag 
(K,< 11.2 mag) 



2.3. Observations 

The two discovered binary systems, 2M 0429 and 2M 1847, were detected at the VLT 
and Subaru observatories on 2003 February 13 (UT) and 2003 July 10 (UT), respectively. 
A total of 6 dwarfs from our sample were observed during these two runs. Table 1 lists the 
4 low- mass dwarfs observed with no likely physical companion detections between ~ O'/l- 
15". For completeness we also include the 27 single stars observed in Paper I. Table 2 lists 
the observable properties of the 2 new binary systems along with the 3 systems presented 
in Paper 1. Target stars were considered "observed" when a minimum corrected FWHM of 
~ 0'.'15 in H band was achieved. 

Each of the observations were made by dithering over 4 different quadrant positions on 
the infrared camera detector. For all targets we obtained both unsaturated H or images 
(< 10s, "short" images), depending on seeing conditions, and saturated H images (30s, 
"deep" images) to gain sensitivity to potential faint companions. 
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At Subaru we used the Coronagraphic Imager with AO (CIAO) without using the 
coronagraphic mask feature. The detector is a 1024x1024 ALADDIN II InSb infrared hybrid 
array with a platescale of 0'.'02 17 pixel" ^ (Tamura et al. 2000). For 2M 1847 we took a total 
of 12xl0s short exposures at H and K^, 12x20 s at J, and 12x60 s deep exposures at H. At 
the VLT wc observed with the Nasmyth AO System/NIR Imager and Spectrograph (NACO) 
system on UT4 (Yepun) which contains a 1024 x 1024 ALADDIN II InSb infrared hybrid array 
detector with a platescale of 0'/0271pixel~^. NACO is unique in that it utilizes an infrared 
WFS. We found that the infrared WFS was most efficient for objects with Kg < 11. 2 mag. 
For 2M0429 we took a total of 16x0. 5s short exposures at H, 8x0. 5s at K^, 12x1 s at J, 
and 12 X 30 s deep H frames. 



2.4. Reduction 

The images were reduced using an AO data reduction pipehne written in the IRAF 
language as first described in Close et al. (2002a). The pipehne produces final unsaturated 
exposures in J, H, and Kg with deep 720 s exposures at H band for each observed binary 
system. The dithering of the shorter exposures produces a final 30" x 30" image with a high 
S/N region in a 10" x 10" box centered on the binary. In order to detect close companions 
within 1" of the central star we filter out the low spatial frequency components of the deep 
images leaving behind high frequency residuals in the PSF (unsharp masking). No faint 
companions, however, were found within the halo of our central stars using this technique. 
Both binary systems were detected from reductions of the shorter exposures. Figures 1 and 
2 show Kg images of the two new systems. 

Photometry for the more widely separated 2M 0429 was performed using the DAOPHOT 
PSF fitting photometry package in IRAF. The PSFs used were unsaturated single stars 
observed during the same night with similar IR brightness, spectral type, and air mass. The 
errors in Amag, listed in Table 2, are the differences in the photometry between 2 similar 
PSF stars. 

DAOPHOT could not successfully separate the strongly blended components of 2M 1847 
due to lower Strehl ratios caused by observing through a 1.4 airmass (the Strehl ratio and 
hence resolution were better at airmass of 1 when the binary was initially discovered and its 
components more clearly separated, however, technical difficulties resulted in delayed image 
acquisition). Consequently, we remove the low spatial frequencies of the binary revealing 
their high-frequency cores. We then perform aperture photometry using IRAF PHOT. This 
purely differential technique preserves the relative magnitude difference between each com- 
ponent while removing sufficient primary halo fiux to reveal the companion centroids in K^. 
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The technique gave reliable Amags and was verified on binary images with known Amags. 

We calculate individual fluxes and their uncertainties from the measured binary flux 
ratios and the integrated 2MASS apparent magnitudes (2MASS All-Sky Point Source Cat- 
alog), along with their respective uncertainties. Table 3 lists the photometry and derived 
characteristics of the new binary systems. 

2.5. An Example of Sensitivity: The Special Case of 2MASS 0253 

One of our targets observed to have no stellar companions deserves special mention. 
2M02530084+1652532 (hereafter 2M 0253) is a newly discovered M6.5 dwarf (Teegarden et 
al. 2003) remarkably only 3.6 pc away. It was discovered from a search of the SkyMorph 
database of the Near Earth Asteroid Tracking (NEAT) project (Pravdo et al. 1999) as a 
high proper motion object (5"yr~^). The star's proximity presented a rare observational 
window for the direct imaging of several Jupiter-mass, extrasolar planets. We were able 
to both probe semimajor axis separations to within ~3AU of the star, comparable to the 
separations of known extrasolar planets detected through radial velocity studies (a< 6 AU; 
http://exoplanets.org) and outside the speckle-dominated region on the detector (> 1"). 
We were the flrst to observe this object with high resolution on 2003 July 14 (UT) using 
the NIRC2 camera and AO (Wizinowich et al. 2000) on the W. M. Keck II telescope. The 
OVOlpixel""^ plate scale mode was used on the 1024x1024 pixel array. 

2M 0253 was only observable for approximately 1 hour in the early morning. We achieved 
sensitivity to companions of 11=19.6 at 1'.'5 in 24min of total integration time (49 x 30 s frames 
in a 4-dither pattern) . We fully saturated the central star so as to allow for the detection of 
any massive faint Jupiter planets orbiting > 1"(> 2.6 AU) from the central star. No faint 
companions were detected. 

The 24min of total integration time enables us to establish upper limits on planetary 
masses orbiting this star. We construct an unsaturated PSF of the star by replacing the 
saturated core with scaled unsaturated pixels from a short exposure. We determine maximum 
H band A mag contrasts by combining scaled models of faint companions (with appropriate 
PSFs) at various radial distances until a 5a detection is obtained. Figure 3 shows the 
resulting 5o" limiting magnitudes at several radial distances from the central star. The 
horizontal dashed lines indicate the H band A mag required for the detection of 5 Gyr, 10 Mj 
and 25 M J objects using the models of Burrows, Sudarsky, & Lunine (2003). We use the 
peak of their H band spectra to estimate the flux emission in this exercise. The star's age 
is not known but based on its high tangential velocity we can assume it is an older object 
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(Wielen 1974). The figure demonstrates sensitivity to an 11 Jupiter-mass extrasolar planet 
at only ~4AU (1.5") away. In Figure 4 we show the fully reduced Keck image of 2M0253 
spatially filtered of its low frequency halo, leaving behind high frequency residuals in the 
core (supcrspcckles). This image also illustrates that with conventional AO, speckle noise 
limits the detection of faint companions within the inner ~ 1" of the halo. 

3. Analysis 

3.1. Are the Companions Physically Related to the Primciries? 

From the total of 69 objects already observed in both this survey and a companion 
survey of M8.0-L0.5 stars by the authors (Close et al. 2003), we did not detect any additional 
unknown red (J-K^ > 0.8 mag) background objects in 6.2 x 10^ square arcsec. Therefore, 
we estimate the probability of a chance projection of a comparably red object within 0.5" 
of the primary to be < 1.3 x 10~^. As we argued in Paper I, with an M6-M8 dwarf density 
of 0.007 pc~^ in the local solar neighborhood (Reid & Gizis 1997a), the probability of an 
apparent companion being just a background star at, for example, twice the distance of 
the target star (hence fainter by a A magnitude of 1.5 mag) and appearing within 0.5" of 
any of our targets is estimated to only be ~ 3 x 10"''. Additionally, none of the companion 
images appear spatially extended as might be expected of background galaxies. Therefore, we 
conclude that both of the very red companions are physically associated with their primaries 
and hereafter we will refer to them as 2M 0429B and 2M 1847B. 



3.2. Distances 

Neither of the 2 binary systems have published trigonometric parallaxes. Wc estimate 
distances to both primaries from a color-magnitude diagram developed in Paper I based on 
trigonometric parallaxes of other well-studied, late-M, field dwarfs from Dahn et al. (2002). 
Using corresponding 2MASS photometry for each star with a trigonometric parallax, we 
estimated a linear least-squares fit of M^^ = 7.65 -|- 2.13 (J-K^) for the spectral range M6.5- 
Ll. This relationship has a la error of 0.33 mag, which has been added in quadrature to 
the J and K^. photometric errors to yield the primarys' Mk^ values reported in Table 3. We 
then use the distance modulus of the primary to estimate the distances to the binaries. The 
calculated distances are listed in Table 3. 
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3.3. Spectral Types and Temperatures 

We do not have spatially resolved spectra of the individual components in either of the 
2 new systems. We estimate the spectral types of each of the binary components by using 
the relation SpT = 3.54M/^^ - 27.20 derived in Paper I from the data set of Dahn et al. 
(2002) (eg. SpT = 8 is an M8, SpT = 10 is an LO, etc). This relationship has a 1 cr error of 
0.85 spectral types which when taken in quadrature with the uncertainty in M^^^ gives an 
overall uncertainty of about 1.5 spectral types. Fortunately, none of analysis is dependent 
on these spectral type estimates. The results are listed in Table 3. 

Effective temperatures of the binary components are estimated from the DUSTY evo- 
lutionary tracks (Chabrier et al. 2000) using calculated M^^ values and estimated ages 
(see Figures 5 and 6). We estimate 2M0429A and 2M0429B to have effective tempera- 
tures of 2690+170 K and 2240l26o K, respectively; 2M 1847A and 2M 1847B are estimated at 
2760^260 ^^'^ 2690+210 -^j respectively. These estimated temperatures are in very good 
agreement with those predicted in Dahn et al. (2002) for the given spectral types. 

3.4. Ages and Masses 

Estimating the age of late-type field dwarfs without Li measurements or established 
cluster membership is difficult. Consequently, we conservatively assume a mean age of ~ 
5 Gyr for our objects with uncertainty spanning the range of common ages in the solar 
neighborhood (0.6-7.5 Gyr; Caloi et al. 1999). 

To estimate masses of these objects we rely on luminosity-mass- age models for VLM 
stars and brown dwarfs. We utilize the DUSTY models to provide theoretical estimates for 
both stellar and substellar masses as a function of both absolute magnitude and age. 
The tracks are calibrated for the K., bandpass (I. Baraffe, private communication) and we 
extrapolate the isochrones from 0.10 to 0.11 M© so as to enclose the upper mass limits of 
our central stars. The companion's absolute magnitude is simply determined by adding the 
measured AK^ to its primary star's M^^- The crosses in Figures 5 and 6 indicate the best 
estimates of where the binary components lie on the 5 Gyr tracks and their uncertainties 
are represented by the shaded regions. 2M 0429A's region of uncertainty is displayed in the 
upper right while its companion is displayed in the lower left. Because the 2M1847 binary 
system is of near equal magnitudes, their regions of uncertainty largely overlap. In this case 
the slightly more massive primary's region of uncertainty is indicated in bold outline and the 
portion of the companion's not overlapping is dashed. The maximum mass is related to the 
minimum M;^^ at the oldest possible age; the minimum mass is related to the maximum Mk^ 
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at the youngest possible age. Table 3 lists the estimated masses for both binary systems. 
Both systems' primary masses arc consistent with M7-type dwarfs and their secondaries 
are most likely stellar, however, 2M0429B's uncertainties extend well into the substellar 
region according to the model. The uncertainty in the masses, as well as in the effective 
temperatures, is largely driven by the uncertainty in our determination of M;^^ (a = 0.33 
mag) as obtained from the [M^^ , J-K^] color magnitude diagram linear fit from Paper I. 
Future observations of trigonometric parallaxes would significantly reduce the uncertainty 
in and hence the masses and temperatures. 



4. Discussion 

4.1. The Binary Frequency of M6.0-M7.5 Stars 

We update the binary fraction statistics of M6.0-M7.5 stars combining the latest results 
presented here (2 binaries resolved out of 6) with those from Paper I (3 binaries out of 30). 
This implies an observed, uncorrected binary fraction of 14l4% using a Poisson distribution 
for the uncertainty (Burgasser et al. 2003). However, this sample was originally drawn from a 
magnitude-limited sample and hence the observed binary fraction is biased due to the leakage 
of equal magnitude binaries into our sample from further distances (Malmquist bias). We 
therefore need to correct our result due to this bias as well as consider sample incompleteness 
due to undetected very tight lower mass companions. 

To compensate for the fainter single stars not inchidcd in our fiux-limitcd sample we first 
adjust for a larger observed volume due to the discovered binaries by a volume correction 
factor. This factor is simply the ratio of the spherical volume containing 95% of our detected 
binaries and the spherical volume containing 95% of our target objects. This results in a 
volume correction factor of (30pc/24pc)^ = 2 and a Malmquist-corrected binary fraction of 
5/(36 X 2) or 7tt%. 

The possibihty that there were faint companions, both stellar and non-stellar, not de- 
tected due to instrument insensitivity is a real one. The curve in Figure 7 shows the in- 
strumental sensitivity of our sample in the speckle noise hmited region (< 1", 30 AU for a 
star assumed 30 pc away) as a function of mass ratio and projected separation in AU. It is 
based on modelling of a 5 Gyr (typical of the ages expected for field stars Gizis ct al. (2000)), 
M6.5 dwarf placed at 30 pc (typical of the distances of our discovered binaries) observed at 
the 8-m Gemini North telescope. We use the DUSTY models to convert AH magnitudes to 
mass ratios. The reason we convert to mass ratios is because it allows us to use the observed 
mass ratio distribution for VLM binaries (Close et al. 2003) to predict the number of missed 
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companions with q> 0.6. The 5 large asterisks in Figure 7 represent the 5 discovered binary 
systems from this survey. Interestingly, they are all found in the upper left corner of the 
sensitivity curve. The fact that some are so near the curve strongly infers that binaries just 
below the sensitivity curve were most likely missed. 

To apply an instrument-sensitivity correction we need to estimate how many binaries 
went undetected in our sample. We generate a Monte Carlo simulation of 11670 synthetic 
companions with the binary properties of VLM systems. For our model we use the mass 
ratio and separation distributions for VLM binaries and assume that the distributions are 
independent. For the mass ratio distribution wc assume a power law decline from unity to 
0.6 from Close et al. (2003). We create the separation distribution profile by plotting the 
42 most tightly separated and resolved VLM {Mtot < 0.19 Mq) binaries currently known 
(see Figure 8; Table 4). Originally presented in Close et al. (2003), we update the list of 
all known VLM binaries and present it in Table 4. The definition of VLM binary having 
a total mass of <O.19M0 is selected to ensure that the binary components are of spectral 
type M6.0 or later and hence differentiated from more massive systems. The peak of this 
distribution, ~ 5AU, is much tighter than the ~ 30 AU distribution peak of slightly more 
massive M0-M4 dwarfs (Fischer & Marcy 1992) and solar-mass stars (Duquennoy & Mayor 
1991). The separation distribution is bound by the smallest separation the instruments were 
able to obtain in H band (~ 0'/08 x 30 pc) on the near side and the empirically sampled wider 
separation encompassing 95% of known VLM binaries (Table 4) on the far side. 

From this sample of nearly 12,000 simulated companions, 21% were below the instrument 
sensitivity curve (but above the instrument sensitivity mass ratio cutoff of q = 0.6) as shown 
in Figure 7. With 5 detected binaries, this predicts 1.3 companions were missed in our 
survey. Hence we conclude that the binary fraction for M6.0-M7.5 stars is (5-|-1.3)/36/2 or 
91*13%. It should be pointed out, however, that the true fraction is most certainly larger than 
this figure since we cannot rule out the possibility of the existence of low q binaries due to 
the sensitivity of this survey. Therefore, our reported binary fraction, accurate within the 
uncertainties for M6.0-M7.5 dwarfs for separations > 3 AU, represents a low-end estimate to 
the intrinsic binary fraction. 

With slightly improved statistics, this latest result for the binary fraction of M6.0-M7.5 
stars is now more comparable with those of later spectral types: 15±7% for late M/early 
L dwarfs (Close et al. 2003), 10-15% and 15±5% for L dwarfs (Bouy et al. 2003; Gizis et 
al. 2003, respectively) and 9li\^% for T dwarfs (Burgasser et al. 2003). These cooler dwarfs 
including the ones presented here all have binary fractions significantly lower than the ~ 32% 
observed for earher M0-M4 dwarfs (Fischer & Marcy 1992) and the ~ 50% for solar-mass 
stars (Duquennoy & Mayor 1991) over the same a>3AU separation range. Our conclusion 
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from Paper I is strengthened: for spectral type M6.0-M7.5 binary systems with separations 

SAU< a< 300 AU the binary fraction from our survey is fff^Yo, statistically consistent with 
cooler M, L, and T stars and significantly less common than that of G and early M stars. 

4.2. The Separation Distribution Function 

The 2M0429 and 2M1847 binary systems have best-estimate projected separations of 
6 AU and 2 AU, respectively. In our total M6.0-M7.5 sample, we detect no binary separations 
wider than 10 AU (Table 3). When we analyze the semimajor axis separations of these 
binaries we observe that there are no projected VLM binary separations > 15 AU. With our 
survey sensitive out to ~ 15" from the central star, this result appears to be real and not 
a result of a sensitivity selection effect. When examining the entire VLM sample of known 
binary systems from the literature (44, Table 4), only three objects are currently known to 
have a projected separation > 15 AU. This indicates that while wide VLM binaries of q > 0.6 
can exist, they are rare. 

The median projected separation of our entire binary sample is ~ 5AU, consistent with 
the ~4 AU peak distribution of late M/early L binaries (Close et al. 2003), L dwarfs (Gizis et 
al. 2003; Bouy et al. 2003), and T dwarfs (Burgasser et al. 2003). This contrasts significantly 
with the ~ 30 AU broad separation peak of early M and G dwarfs (Fischer & Marcy 1992; 
Duquennoy & Mayor 1991). We conclude that the projected semimajor axes of M6.0-M7.5 
binaries appear consistent with those of late M, L, and T dwarf systems but are significantly 
smaller on average than early M and G stars. 

5. Summary 

We have conducted the largest flux hmited (K^ < 12 mag) survey of nearby M6.0-M7.5 
dwarfs using the Keck II, Gemini North, Subaru, and VLT AO systems. In this paper 
we present our 2 latest binary discoveries, 2M 0479 and 2M 1847, observed at the VLT 
and Subaru facilities, respectively. When added to our initial results from Paper 1, the 
overall survey consists of 36 stars with 5 discovered binary systems. The 2 new components 
are of relatively equal mass (q>0.8) with average projected separations of 2 and 6AU. 
While none of the binaries have yet been confirmed by common proper motions, they are 
almost certainly bound based on space density arguments of very red companions. We have 
used observational and statistical arguments to characterize the VLM binary fraction and 
separations that contribute additional empirical constraints to binary formation mechanisms: 
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• We estimate the binary frequency of spectral type M6.0-M7.5 main sequence stars for 
separations a>3 AU from this survey to be 9^^%. The figure is statistically consistent 
with later type ultracool M, L, and T dwarfs. The frequency is significantly less than 
that measured in studies of earlier M and G dwarfs, inferring that the binary fraction 
of stars is a function of the spectral type of the central star. 

• The separations of the 5 binary systems from our sample are all < 10 AU. Projected 

separations of known VLM binaries > 15 AU arc rare. This survey's median separation 
of 5 AU is consistent with the separations of later type M, L, and T dwarfs (separation 
peak ~ 4AU). This is in stark contrast with the broad peak separations of ~ 30 AU 
for the more massive M and G binaries. 
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Table 1. M6.0-M7.5 Stars Observed with No Likely Physical Companion Detections 

Between 0.1"-15"^ 



2MASS Name 


Other Name 


Ks 


Spectral Type 


Reference 


2MASS J0253008+165253'' 


SO 025300.5+165258 


7.59 


M6.5 


5 


2MASSI J0330050+240528 


LP 356-770 


11.36 


M7.0 


1 


2MASSI J0435161-160657'' 


LP775-31 


9.34 


M7.0 


2 


2MASSI J0752239+161215 




9.82 


M7.0 


2 


2MASSI J0818580+233352 




11.13 


M7.0 




2MASSW J0952219-192431 




10.85 


M7.0 




2MASSW J1016347+275150 


LHS 2243 


10.95 


M7.5 




2MASSI J1024099+181553 




11.21 


M7.0 




2MASSW J1049414+253852 




11.39 


M6.0 




2MASSI J1124532+132253 




10.03 


M6.5 




2MASSW J1200329+204851 




11.82 


M7.0 




2MASSW J1237270-211748 




11.64 


M6.0 


J 


2MASSW J1246517+314811 


LHS 2632 


11.23 


M6.5 




2MASSI J1253124+403404 




11.20 


M7.5 




2MASSW J1336504+475131 




11.63 


M7.0 




2MASSW J1344582+771551 




11.83 


M7.0 




2MASSI ,11356414+434258 




10.63 


M7.5 


2 


2MASSI J1431304+171758'' 




11.16 


M6.5 


2 


2MASSI J1521010+505323'' 




10.92 


M7.5 


2 


2MASSP J1524248+292535 




10.15 


M7.5 


3 


2MASSW J1527194+413047 




11.47 


M7.5 


3 


2MASSW J1543581+320642 


LP 328-36 


11.73 


M7.5 


1 


2MASSW J1546054+374946 




11.42 


M7.5 


1 


2MASSW J1550381+304103 




11.92 


M7.5 


1 


2MASSW ,11757154+704201 


LP 44-162 


10.37 


M7.5 


1 


2MASSW ,12052086-231809 


LP 872-22 


11.26 


M6.5 


1 


2MASSW J2221544+272907 




11.52 


M6.0 


1 


2MASSW J2233478+354747 


LP 288-31 


11.88 


M6.0 


1 


2MASSI J2235490+184029 


LP 460-44 


11.33 


M7.0 


1 


2MASSW J2306292-050227 




10.29 


M7.5 


1 


2MASSW J2313472+211729 


LP 461-11 


10.42 


M6.0 


1 



Tor near-equal mass companions. For smaller companion masses with q <0.8, sensitivity is a 
function of distance. See Figures 3 and 7. 

''Results from this paper; otherwise, Paper L 



References. — (1) Gizis et al. (2000) (2) Cruz et al. (2003) (3) Reid et al. (2002) (4) Kirkpatrick, 
Henry, & McCarthy (1991) (5) Teegarden et al. (2003). 
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Table 2. The New Binary Systems 



System 


AJ 


AH 


AXa 


Sep. (mas) 


P.A. (deg) 


Date Observed (UT) 


Telescope 


LP 415-20" 


0.84 ± 0.15 


0.77 ± 0.10 


0.66 ± 0.06 


119 ± 8 


91.2 ± 0.7 


2002 Feb. 07 


Gemini North 


LP 475-855*' 


0.48 ± 0.05 


0.43 ± 0.04 


0.48 ± 0.03 


294 ± 5 


131.6 ± 0.5 


2001 Sep. 22 


Gemini North 


2MASSI .10429184-312356= 


1.20 ± 0.12 


1.10 ± 0.08 


0.98 ± 0.08 


531 ± 2 


298.9 ± 0.2 


2003 Feb. 13 


VLT 


2MASSW J1750129-H442404 


0.74 ± 0.15 


0.73 ± 0.15 


0.64 ± 0.10 


158 ± 5 


339.6 ± 0.7 


2002 Apr. 25 


Gemini North 


2MASSI J1847034-|-BB2243'= 


0.26 ± 0.18 


0.34 ± 0.15 


0.16 ± 0.10 


82 ± 5 


91.1±1.4 


2003 July 10 


Subaru 



^Also known as Bryja 262. 

^'Also known as [LHD94] 042614.2+13312 and 2MASSW J0429028+133759. 
'^Results from this paper; otherwise, Paper I. 



Tabic 3. Summary of the Binaries' A and B Components 



Name J H Ks Mjf^ * SpT*' d^y^^t {pc)" Mass (Mq)'' Sep. (AU)" P (yr)*" 



LP 416-20A 


13, 


,09 


± 


0.06 


12, 


,47 


± 


0, 


,05 


12, 


12 


± 


0.04 


9.72 


± 


0.38 


M7.0 


30 


± 


5 


0, 


,097 


+ 0,011 
-0.012 


3.6 


± 


0.7 




LP416-20B 


13. 


,93 


± 


0.16 


13, 


.24 


± 


0, 


11 


12, 


,78 


± 


0.08 


10.37 


± 


0.39 


M9.5 








0, 


,081 


+ 0.009 
-0,010 










LP475-855A 


13. 


21 


± 


0.04 


12, 


.54 


± 


0, 


,04 


12, 


18 


± 


0.04 


9.84 


± 


0.36 


M7.5 


29 


± 


5 


0, 


,093 


+ 0,012 
-0,009 


8.6 


± 


1.6 




LP475-855B 


13. 


,69 


± 


0.07 


12 


,97 


± 


0, 


,06 


12, 


,66 


± 


0.05 


10.32 


± 


0.36 


M9.5 








0, 


,082 


+ 0,009 
-0,009 










2M 0429AS 


11, 


18 


± 


0.04 


10 


55 


± 


0, 


,03 


10, 


14 


± 


0.03 


9.88 


± 


0.35 


M7.5 


11 


± 


2 


0, 


,094 


+ 0,010 
-0,011 


6.0 


± 


1.0 




2M 0429B8 


12, 


,38 


± 


0.13 


11 


,65 


± 


0, 


09 


11, 


12 


± 


0.07 


10.86 ± 0.36 


Ll.O 








0, 


079 


+ 0,005 
-0,018 










2M 1750 A 


13, 


23 


± 


0.06 


12 


62 


± 


0, 


06 


12, 


24 


± 


0.05 


9.77 


± 


0.39 


M7.5 


31 


± 


6 





()97 


+ 0.012 
-0.012 


4.9 


± 


0.9 




2M1750B 


13, 


97 


± 


0.16 


13 


35 


± 


0, 


16 


12, 


88 


± 


0.11 


10.41 


± 


0.41 


M9.5 











085 


+ 0.006 
-0.016 










2M 1847AS 


12, 


55 


± 


0.08 


11 


87 


± 


0, 


07 


11, 


,58 


± 


0.05 


9.72 


± 


0.41 


M7.0 


23 


± 


4 


0, 


,098 


+ 0.022 
-0.012 


1.9 


± 


0.4 


Q+3 


2M 1847B'5 


12, 


,81 


± 


0.20 


12, 


,21 


± 


0, 


16 


11, 


,74 


± 


0.11 


9.88 


± 


0.42 


M7.5 








0, 


,094 


+ 0.014 
-0.013 











^M^^ ^ 7.65 + 2.13(J--fCs) with a rms <y m ^ = 0.33 derived in Paper I (§4.2); relationship is valid for M6.5<SpT<Ll. 

'^Spectral type estimated by SpT — 3.54Mj^, -27.20 with itl.5 speetral subelasses of error in these estimates as derived in paper I (§4.5). SpT = 
10 is defined as an LO; valid for M(i.5<SpT<Ll. 

'^Distances based on M/^ as described in §4.2. 

'^Mass determination uses the models of Chabrier et al. (2000). 

^Projected separations. 

^Periods include a 1.26 multiplication of the projected separations compensating for random inclinations and eccentricities (Fischer &c Marcy 
1992). 

®This paper; otherwise, Paper I. 
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Table 4. All Known Resolved VLM Binaries'" 



Name 


Sep.'' 


Est. 


Est. Ma 


Est. Mb 


Est. Period'' 


Ref.' 




AU 


SpTa/SpTb 


Mq 


Mq 


yr 




PPL 15 


0.03 


TV T ^ / A T O 

M7/M8 


0.07 


0.06 


5.8 days 


1 


Gl 5d9B 


1.0 


TV JO /TV X r\ r\ 

JV18.5/M9.0 


0.063 


0.06 


3 


2,3 


bUbb ioJo5oj-UUloU4 


1.1? 


L17/L4? 


0.079 


0.074 


3 


4 




1.5 


L4/L6 


0.073 


0.070 


5 


4 




1.8 


1 0.5/ 1 0.0 


0.05 


0.05 


8 


5 


2MASSI J1847034+552243 


1.9 


M7/M7.5 


0.098 


0.094 


9 


This paper 


iMAsbW J085D479+2i3518 


2.0 


T r O It aO 

Loi /LSc 


0.071 


0.064 


8 


4 


Dii.iNib-F J 185950. y-ovUboz 


2.0 


L0/L3 


0.084 


0.076 


7 


4 


riU loUy4ors 


2.4 


L2/L2 


0.07 


0.06 


10 


6 


e incli o 


2.0 


1 1/ lb 


U.U42 


U.U27 


Id 


1 

lb 


OUT A OOTT17 TO^^/^^or' 1 OfAOOOO 

2MAbbW J074d425+200032 


2.7 


L0/L1.5 


0.085 


0.066 


11 


7, 17,20 


zMAobW J1U47127+4U2d44 


2.7 


JVL8/1jU 


0.092 


0.084 


11 


8, 17 


DrL(JNiD-F J(Jo572b.9-4417oU 


2.8 


L2/L4 


0.078 


0.074 


12 


4,13 


oA/TAooTsr Tioocrc/iQo orTin/i^^r* 

2MAboW J i225o4o2-27o94bD 


3.2 


1 b/ 1 5 


0.033 


0.024 


23 


5 


zMAbo W J 0921)122+051742 


3.2 


L6.5/L7 


0.068 


0.068 


16 


7 




0.4 


A /TQ c; /TV TO n 
M8.0/iVi9.U 


0.070 


0.068 


16 


18 


415-2U 


3.5 


TV Tf7 /Tl TO r- 

M7/M9.5 


0.095 


0.079 


15 


9, Paper I 


on T A OOTX7 T1 TOO 1 t A l QO/lOCn 

2JVlAbbW J172oll44-o94o59 


3.7 


L7/L8 


0.069 


0.066 


19 


4,13 


T TTO OOOT 

LHb 2397a 


3.9 


TV TO /T T r" 

M8/L7.5 


0.090 


0.068 


22 


10,17 


zMAboW J 14zoolD+loO /Ul 


3.9 


Mo.O/ IjI 


0.088 


0.076 


19 


17 


2JVlAboW J214U29o+lo251o 


3.9 


A Tn /T o 

M9/L2 


0.092 


0.078 


22 


17 


zMAooW JzZUozZo-zU4 ( Uo 


4.4 


A (T Q /A /TO 

Mo / Mo 


0.092 


0.092 


22 


1 7 

17 


2JVlAbbs JUooUo59+1Uo71d 


4.4 


L6/L8 


0.05 


0.04 


30 


7 


2MAboW J175U129H-4424U4 


4.8 


A /TT cr /T n 

M7.5/LU 


0.095 


0.084 


25 


9, Paper i 


■pvTTiTVTTO T~) T1000 O 1 C AT 

Dhj JN lb- r* J 1228.2-1547 


4.9 


L5/L5 


0.05 


0.05 


34 


11 


OA T A 0011(7 T1 cr\nr\cr a i i taooo 

2MAbbW J1dUUU54+17U8o2 


5.0 


L1/L3 


0.078 


0.075 


29 


4,13 


zMAsbW J12j9272+5515o7 


5.1 


L5/L5 


0.071 


0.071 


31 


4,13 


zMAbbl JU429184-jli35D 


6.0 


TV /T 'T r /T 1 

M7.5/L1 


0.094 


0.079 


50 


This paper 


IFJVIBD 29 


7.2 


Ll/L4^ 


0.045 


0.038 


68 


14 


zivi/\oo vv J 1 i^Do^o-rZZouoo 


1 .o 


T /T A 


u.uoo 


u.uoo 


Do 


1 9 


2MASSW J131 1391+803222 


7.7 


M8.5/M9 


0.089 


0.087 


51 


17 


CFHT-Pl-12 


7.8 


M8/L4S 


0.054 


0.038 


76 


14 


2MASSWJ1127534+741107 


8.3 


M8/M9 


0.092 


0.087 


57 


17 


LP 475-855 


8.3 


M7.5/M9.5 


0.091 


0.080 


58 


9, Paper I 


DENIS-PJ0205.4-1159 


9.2 


L7/L7 


0.07 


0.07 


75 


15 


2MASSW J2101349+175611 


9.6 


L7/L8 


0.068 


0.065 


82 


4,13 


2MASSW J2147436+ 143131 


10.4 


L0/L2 


0.084 


0.078 


83 


4,13 


2MASSW J1449378+235537 


11.7 


L0/L3 


0.084 


0.075 


100 


4,13 


IPMBD 25 


11.8 


M7/L4g 


0.063 


0.039 


126 


11 


DENIS-P J144137.3-094559 


13.5 


Ll/Ll 


0.079 


0.079 


124 


4,13 


2MASSWJ2331016-040618 


15.0 


M8.5/L7 


0.093 


0.067 


159 


17 


CFHT-Pl-18 


34.5 


M8/M8 


0.09 


0.09 


641 


4,19 


2MASSWJ1207334-393254 ^ 


55 


M8/L5-L9.5 


0.02 


0.005 


2600 


22 


2MASSW J11011926- 7732383 


240 


M7/M8 


0.05 


0.025 


13600 


21 



''We define VLM binaries as systems whose total mass is < 0.19 M©. Very young evolving systems like GG 
TauBaBb (White & Ghez 2001) are not included, nor are overluminous systems that are not resolved into 
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binaries. 

''Projected separation except for the systems where semimajor axes have been measured. 

•^This "period" is simply an estimate assuming a face-on circular orbit except for the systems where P has 
been measured. 

''PPL 15 is a spectroscopic binary (Basri & Martin (1999) and is not resolved). 
These tight binaries have a widely separated more massive primary. 
^Part of a triple system. 

s Spectral type estimated from the given masses using Pleiades age, effective temperatures from Burrows et 
al. (2001), and spectral type-temperature relation for ultracool dwarfs from Dahn et al. (2002). 

''Parallaxes and common proper motions have yet to be confirmed 

'REFERENCES-(l)Basri & Martin (1999); (2) Kenworthy et al. (2001); (3) Lane et al. (2001); (4) Bouy et 
al. (2003); (5) Burrows, Sudarsky, & Lunine (2003); (6) Potter et al. (2002); (7) Reid et al. (2001); (8) Reid et 
al. (2002); (9) Siegler et al. (2003); (10) Freed, Close, & Siegler (2003); (11) Martin et al. (1999); (12) Koerner 
et al. (1999); (13) Gizis et al. (2003); (14) Martin et al. (2003); (15) Delfosse et al. (1997); (16) McCaughrean et 
al. (2004); (17) Close et al. (2003); (18) Leinert et al. (2001); (19) Martin et al. (2000); (20) Bouy et al. (2004); 
(21) Luhman (2004); (22) Chauvin et al. (2004) 
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Fig. 1. — An 8 X 0.5 s image of the newly discovered binary system 2M0429 shown in the 
K, band; observed on 2003 February 13 (UT) at the VLT. The platescale is 0". 02 71 pixel" ^ 
The contours are linear at the 80, 60, 40, and 20% levels; north is up and east is to the left. 
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a) 2M1 B47 


b) PSF star 






HC > 


J- - "~ "v. 






_ ^ 





Fig. 2. — (a) A 12 X 10 s image of the newly discovered binary system 2M 1847 shown in the 
band; observed on 2003 July 10 (UT) at Subaru. The platescale is 0".0217pixer^ The 
contours are hnear at the 80, 60, 40, and 20% levels, (b) For comparison, the PSF star 
2M 0253 observed in the same evening and at the same airmass is displayed. The contours 
are linear at the 85, 75, 65, 55, 45, 30, and 15% levels. In both images, north is up and east 
is to the left. 
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Fig. 3. — Instrument sensitivity curve showing 5a AH detection versus distance in arcsec 
from the very nearby M6.5 star 2M 0253 . Total integration time is 24 min using Keck II AO. 
The "crosses" indicate the 5a sensitivity hmits of our data to simulated faint companions. 
The upper horizontal dashed line corresponds to a 5 Gyr old, 25 Mj brown dwarf using 
the models of Burrows, Sudarsky, & Luninc (2003) while the lower horizontal dashed line 
corresponds to a 5 Gyr old, 10 Mj planet. We were sensitive to the detection of a 10 Mj 
planet at 4AU and a 25 Mj brown dwarf at 1.5 AU (assuming 2M0253 is only 3.6 pc away 
and 5 Gyr old) . 
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Fig. 4.— Reduced 24min image of 2M0253 observed at Keck II on 2003 July 14 (UT) 
with its low spatial frequencies removed (unsharp masked). The simulated 5(7 companion 
is 33,000 times fainter than the central star (AH=11.3) at only 3.6 AU (1-5"). Teegarden 
et al. (2003) report H=7.9 for the central star which allows detection to a 5 Gyr old, 11 Mj 
planet according to the models of Burrows, Sudarsky, & Lunine (2003). Also visible are the 
residuals from the 6 spider arms and super speckles. To within our sensitivity limits we find 
no companions to 2M 0253. 
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Fig. 5. — Chabrier et al. (2000) DUSTY stellar and substellar evolutionary tracks custom 
integrated over the bandpass ([m/H] = 0). The best-guess values of the individual binary 
components of 2M 0429 are indicated by the bold "crosses" with the primary at the top right 
and the companion lower and to the left. The shaded polygons enclose each components' 
region of uncertainty. The components' derived M^^ is listed in Table 3. With no knowledge 
of the binary's age, we conservatively assign a mean age of 5 Gyr and uncertainties spanning 
the range of ages in the solar neighborhood (0.6-7.5 Gyr; Galoi et al. 1999). The model 
suggests a primary mass of 0.094^q q}° ^ temperature of 2690^}yo K. For the com- 

panion, the model predicts a mass of 0.079^g o58 ^0 ^ temperature of 22401^260 The 
isochrones plotted are (left to right) 0.6, 0.65, 0.7, 0.85, 1.2, 1.7, 3.0, 5.0, 7.5, and 10.0 Gyr 
(the oldest 4 isochrones are indistinguishable at the given scaling). 
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Fig. 6. — As in Figure 5, but for 2M 1847. In this case the shaded regions of uncertainty of 
the 2 components overlap. Hence, we outhne the primary star's in sohd and the companion 
in dotted hnes (where the two regions overlap solid takes preference). The model suggests 
a primary mass of 0.0981o;q^2 and a temperature of 276OI260 ^'^^ the secondary the 
model suggests a mass of 0.094lo:oi3 M© and temperature of 2690l2?S K- 
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Fig. 7. — The results from a Monte Carlo simulation generating 11,670 companions dis- 
tributed according to a bivariate distribution (mass ratio q and separation a) is plotted over 
the instrumentation sensitivity curve (connected lines). We assume the two distributions 
are independent. We assume a power law declining from unity to 0.6 for the mass ratio 
distribution (Close et al. 2003) and the profile from Figure 8 for the separation distribution 
for a > 3 AU. The instrumentation sensitivity curve is based on modeling of a ~ 5 Gyr M6.5 
dwarf placed at 30 pc, typical of the distances of our discovered binaries. The DUSTY models 
(Chabrier et al. 2000) are used to convert AH magnitudes to mass ratios. The 5 discovered 
binary systems are indicated by large asterisks. 21% of the synthetic companions fall below 
the instrumentation sensitivity curve but above the instrument sensitivity mass ratio cutoff 
of q=0.6. This results in a sensitivity correction of 1.3 binaries. 
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Fig. 8. — A histogram of 42 VLM {Mtot < 0.19 Mq, component spectral types >M6.0) 
binaries from Table 4 are plotted (for reasons of clarity, we leave out the 2 widest systems) . 
The distribution is incomplete less than ~ 3 AU; Poisson error bars are plotted. The declining 
profile is a real feature of the distribution as is the paucity of wide binaries greater than 
15 AU. These features along with a tighter peak distribution (~ 3 — 4AU) are significantly 
different from those of more massive primary stars (Fischer & Marcy 1992; Duquennoy & 
Mayor 1991). 



